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We investigate the dynamics of a single, overdamped colloidal particle, which is driven by a 
constant force through a one-dimensional periodic potential. We focus on systems with large barrier 
heights where the lowest-order cumulants of the density field, that is, average position and the 
mean-squared displacement, show nontrivial (non-diffusive) short-time behavior characterized by 
the appearance of plateaus. We demonstrate that this "cage-like" dynamics can be well described 
by a discretized master equation model involving two states (related to two positions) within each 
potential valley. Non-trivial predictions of our approach include analytic expressions for the plateau 
heights and an estimate of the "de-caging time" obtained from the study of deviations from Gaussian 
behaviour. The simplicity of our approach means that it offers a minimal model to describe the 
short-time behavior of systems with hindered dynamics. 

PACS numbers: 05.40.-a, 05.60.-k, 02.50.Ey 



I. INTRODUCTION 

The transport of particles in modulated potentials 
plays a fundamental role in diverse fields such as colloidal 
particles on topologically or energetically structured sur- 
faces GlUf, particles in optical lattices 3| and optical line 
traps [4j , biased Josephson junctions [5j , and in biophysi- 
cal processes [1,0] • In many of these cases, the dynamics 
can be described as overdamped (i.e., non-inertial) Brow- 
nian motion in one spatial direction. A generic model in 
this context is a particle in a periodic "washboard" (e.g., 
sinusoidal) potential, which is tilted by the constant bi- 
asing force, see Fig. [TJt. The resulting particle current 
in the long-time limit was calculated analytically more 
than 50 years ago by Stratonovich || HI • More recently, 
Reimann et al. [lfj derived an analytical solution for 
the corresponding diffusion coefficient related to the long- 
time limit of the mean-squared displacement, (Ax 2 (t)). 
This diffusion coefficient in the modulated potential can 
become much larger than that in a free system, contrary 
to the intuitive idea that confinement slows down diffu- 
sion. This "giant diffusion" effect has also been found in 
experiments [TT| - [l3| . Taken together, one can state that 
the long-time properties of particles in periodic forces are 
nowadays well understood. 

Much less is known on the short-time dynamics of 
driven (or undriven) overdamped particles in periodic 
forces [lj]. However, it is precisely this regime which 
has recently attracted significant theoretical and exper- 
imental interest (see, e.g. 0, )• One reason for 
this interest is the similarity of the observed short-time 
dynamics with that of systems approaching a glass tran- 
sition (l7l-l2l| or, more generally, a state characterized 
by complex energy landscapes yielding trapping. A key 
quantity in this context is the mean-squared displace- 
ment as function of time. For (non-interacting and in- 
teracting) colloids in modulated potentials, this quantity 
has already been investigated by computer simulations 
of the corresponding Langcvin equations (l7| . numeri- 




FIG. 1. (a): The potential landscape of a particle in a tilted- 
washboard. The period is a and the minimum peak-to-trough 
height is Au. (b): The discrete two-state-per-well (2S/W) 
model in which the state of the system is described by well 
index n and a = L,R, describing whether the particle is found 
in the left or right side of the well. The particle therefore 
assumes discrete positions at x — na + x L / R , where x a is 
the displacement of internal state a. Transitions between the 
states are described by a rate equation with rates 7 ± and r 
as shown. 



cal solutions of the Fokker-Planck (Smoluchowski) equa- 
tion [22j, as well as in experiments 0. Provided that 
the barrier height is sufficiently large, (Ax 2 (t)) typically 
grows linearly with time initially, reflecting diffusive mo- 
tion within one valley. At intermediate time a plateau 
appears signaling that the particles are trapped (sub- 
diffusive behavior). Only at long times, do the parti- 
cles eventually overcome the barrier, and (Ax 2 (t)) ap- 
proaches linear time-dependence (these three dynamic 
regimes can be seen in Fig. [3]). These dynamic features 
of a single particle in a modulated potential resemble the 
single-particle dynamics of various many-particle systems 
in kinetically trapped or glassy states [23l 1241 . examples 
being dipolar particles trapped in chains [25j, colloidal 
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gels [26|], and (binary) Lennard- Jones-glasses [23(. In- 
deed, from the perspective of glasses (or gels), the mod- 
ulated one-particle potential is thought to arise from in- 
teraction effects, that is, the formation of cages formed 
by the neighboring particles 17J . Hopping over a barrier 
then corresponds to the escape from the cage. 

Contrary to the long-time properties discussed above 
there is, so far, no simple analytical approach towards the 
transport properties in the short-time regime. Attempts 
have been made to derive, on an approximate level, ex- 
plicit formulae predicting the height of the plateau, as 
well as the crossover time from the sub-diffusive into the 
diffusive regime 0,01 . This yields good estimates, but is 
clearly far away from detailed, quantitative predictions. 
To overcome this gap, we propose in the present paper a 
discretized model, which we treat by the master equa- 
tion formalism. Clearly, the investigation of diffusion 
problems via discretized models and master equations 
has a long tradition (27}. In particular, it is well known 
that conventional master equation approaches without 
memory effects yield a linear time dependence of all cu- 
mulants, i.e., normal diffusive behavior. More complex 
dynamics, also resembling that of glassy systems, is dis- 
played by models based on continuous time random walks 
IialMlilQ and fractional Brownian motion pljlp] 

Here we introduce a discrete model in which the par- 
ticle can adopt one of two positions ("states") within 
each valley. Transitions both between these states and 
between wells are described by rates. Using techniques 
similar to those used in e.g. the theory of full count- 
ing statistics in quantum transport [32H34| . we calculate 
the cumulant generating function for the particles posi- 
tion within this rate equation approach. The treatment 
is simple, requiring just manipulation of two-by-two ma- 
trices, and analytic results can be obtained. In parallel, 
we calculate the density distribution and the cumulants 
numerically from the Smoluchowski equation (SE) and 
demonstrate that our "two-state-per-well" (2S / W) model 
gives an excellent account of the short-time dynamics of 
interest here. In particular, we obtain analytic expres- 
sions for the plateau heights, as well as estimates of the 
"de-caging" or "cross-over" time from the study of the 
nonGaussian parameter, the kurtosis. 

The paper is structured as follows. In section |H] we 
describe the full Smoluchowski equation with a periodic 
potential and in section Hill we present our rate equation 
model. The three dynamical regimes are investigated in 
section II VI and in section [V] we describe the microscopic 
evaluation of the model parameters. Section|VT]compares 
the results between model and numerically exact results. 
Finally, in section fVTTl we conclude. 



II. THE SMOLUCHOWSKI EQUATION WITH 
A PERIODIC POTENTIAL 



We are interested in the diffusion of a particle sub- 
jected to a periodic force in ID. For our purposes the 



full description of the system is given by the probabil- 
ity distribution function p[x, t), which we take to be 
governed by the Smoluchowski (over-damped Fokker- 
Planck) equation Q 



p{x,t) = Cp(x,t), 



with the differential operator 



Do 
kT 



[kTdl + u"{x)+u'{x)d x 



(1) 



(2) 



where kT is the thermal energy, Do the so-called short- 
time diffusion coefficient, and u(x) the modulated poten- 
tial. 

In this work we shall consider as example a tilted wash- 
board potential (our approach is general, however). The 
potential reads 



u(x) = Mo sin ( — ) — Fx, 
a 



(3) 



with a the period of the potential, and parameters uo and 
F describe the depth and tilt of the potential respectively. 
A section of this potential is plotted in Fig. [T^,. 

We assume that the form of the potential is such that it 
may be meaningfully divided up into a number of wells, 
defined between subsequent maxima, and labelled n = 
.... —2, — 1, 0, 1,2,.. .. This occurs when the difference 
between the potential maxima and maxima, Am, is large 
compared with the thermal energy of the particles: Am ^> 
k B T. 

The main quantities of interest here are the cumulants 
of the particles position, which we denote (x k ) c for the 
fcth cumulant. Of particular interest is the second cumu- 
lant or variance (x 2 ) c = (Ax 2 ) = (x 2 ) — (x) 2 . We will 
also consider the kurtosis 



K(t) 



(4) 



as a measure of the deviations from a Gaussian distribu- 
tion. We note that the kurtosis is trivially related to the 
a-parameter often studied in the context of glass forma- 
tion (see, e. g. Ref. [13]): a(t) = §«(*). 



III. TWO-STATE-PER-WELL MODEL 

Our aim is to model the SE for the probability density 
defined as a function of continuous variable x with a sim- 
pler model involving probabilities for finding the particle 
in one of a set of discrete states, defined by regions of the 
x axis. The most obvious approach of this type would be 
to associate one state with each well of the potential and 
write down a rate-equation for the probability (t) of 
finding the particle in well n at time t. This one-state- 
per-well approach is discussed in Appendix [S] and it is 
clear that all position cumulants are linear in time for all 
time. Such a simple rate equation is therefore inadequate 
to describe the short-time dynamics of interest here. 
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The model that we shall focus on is the next most com- 
plicated model in this family, and possesses two states 
per well (2S/W), which correspond to whether the par- 
ticle is to be found in the left or right side of the well, 
see Fig. [TJ). In this model, the location of the particle 
is described by two discrete indices: the integer well in- 

(n) 

dex, n, and the internal index a — L, R and p a is the 
corresponding probability. For calculating the position 
cumulants, we take the particles to be localised at posi- 
tions x = na + x^/r, where x^ and xr are the offsets of 
the left and right states respectively. 

I 



Let us assume that transitions between these states 
can be described by rates and define rates 7 + and 7" 
to describe hopping within a single well, and rates r + 
and r~ to describe hopping between the wells. In both 
cases, superscript + indicates movement to the right and 
— , movement to the left. The rates are homogeneous due 
to the periodicity of the force and, since we assume that 
the barriers are high (Ait 3> kT), the intra- well rates 
7 ± are significantly larger than the inter- well rates T^. 
The system can then be described by the rate equation 
system 






(5) 



Denoting the vector of probabilities with the notation 
n 
3 



p v °')), t his ma ster equation can be written in matrix 



,(«) 
form as [3 



m\p {n) )) + w+\p (n - L> )) + w-\p {n+1> )),(6) 



>-i) 



where Wo,± are two-dimensional matrices defined accord- 
ing to Eq. ([5]). We then define the Fourier-transformed 
vector \p(x;t))) = E n e mx \P {n) (t))), such that, by 
Fourier-transforming Eq. ([6]), we arrive at 

d \p(x;t))) = W( x )\p( X ;t))), (7) 



dV 



with 



W(x) = Wo + e iX W+ + e- lx W-. (8) 
Assuming that the particle starts off in the n = well 
with distribution \p {0) (0))) = {pf\o) , Pr } (0)j , the solu- 
tion of this equation is 

|p(x;t)))-e w ^|p(0;0))). (9) 

Based on this solution, the cumulant generating func- 
tion (CGF) of the particle position is 

F( X ;t) = log((^ | e ^ A e w W t |p(0;0))) (10) 

with internal position matrix 



a 



XL 

x R 



(11) 



and where multiplication from the left with ((0o| = (1, 1) 
traces over internal states [40]. Differentiation of the 
CGF w.r.t. x and taking the x — > limit generates 
the cumulants: 



d(ixr 



(12) 



A proof that the CGF of Eq. (jlQ|) does indeed give these 
cumulants is given in Appcndix[Bj In the following we will 
always start the particle localised within a single well at 
cither position xl or xr. 



IV. LIMITS 

Calculating the CGF just involves manipulation of 
two-by-two matrices and various limiting cases can easily 
be extracted. The three distinct phases of time evolution 
can be looked at in turn. 



A. Short times: linear growth 

Since the intra-well rates 7 ± are larger than the inter- 
well rates T^, in the short-time limit all the dynamics 
occurs within the initial well. We therefore set T ± = 
in the CFG of Eq. (TTUJ) to obtain 



F{x\t) =log< 



b \e ixA e Wi ' 



f |p(0;0))) (13) 



with 



- 7 . 
-, + 



(14) 



To assess the initial movement of the particle, we can 
expand this result to lowest order in t. Of particular 
importance is the evolution of the first two cumulants 
from initially localised distributions at Xljr. These read 



(x(t)) - x L/R ± 7 ± (^i? - x L )t; 



(15) 



where the upper- (lower-) sign corresponds to an initial 
position xl {xr). 
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B. Intermediate times: plateaus 



they should be. The first two read 



Going beyond first order in time, Eq. (|13|) predicts a 
first cumulant, for example, which reads 



(x(t)) c 



7 %L 



l + x R 



7+ + 7 

7 +pi 0) (0)-7-?4 0) (0) 
7+ + 7~ 



{x L - x fl )e-^.(16) 



The first cumulant therefore simply relaxes exponentially 
with rate 7 = 7 + +7~ to an asymptotic value within the 
initial well. The other cumulants behave similarly, with 
the expression for the A:th cumulant being composed of 
a sum of terms with time-dependences e - ™ 7 * where in- 
dex m is an integer with maximum value k. Asymp- 
totically then, the cumulants within a single well reach 
values which correspond to the plateau values. For the 
first four cumulants, the plateau values read 



(x) 



plat _ 7 + ZL+7 XR 



7+ + 7 



fr 2 )g lat = , T^—^ {x L -x R f 
(7+ +7 ) 

(^piat _ 7+7- (7+ -7-) ^ ^ 3 



(17) 



( 7+ + 7 -)3 



-(XL ~ XrY 



4splat _ 7 + 7 (7 ^ 47+7 + 7+ ) 



(x L - X R ) 4 . 



These are the "caged values" . On time-scales larger than 
(r ± ) \ the particle begin to leave the first well and the 
values move away from the above. 



C. Long time-limit 

In the asymptotic t — > 00 limit, the complete CGF of 
Eq. (Til))) becomes 

(Xi t) = Hm JHx; t) = X (x)t (18) 

t— >oo 

where Ao(x) is that eigenvalue of W(x) which behaves as 
lim x ^ Ao(x) = 0: 



Ao(x) 



Tt 

y 



_£ 7 ±r±(e=F*-l) 



, (19) 



with r = 7 + + r + + 7 _ + r~ . From this expression it 
is clear that, in the long-time limit, intra-well distances 
are irrelevant implying that (x k ) c ~ a k (n k ) c . The reason 
is that, in this limit, the extension of the probability 
distribution function is much larger than a single well. 
Moreover, combining Eq. (fT2"j) and Eq. (fT5|). one finds 
that the cumulants are then simply linear in time, as 



(x(t)} 



asym 



7+r+ -7-r- 



(20) 



2 (i)) asym — 



r 2 (7+r+ + 7 _ r~) - 2( 7 +r+ - 7-r- 



(21) 



The latter gives directly the long-time diffusion coeffi- 
cient 

D = lim ~(x 2 (t)) c = ^(x 2 (t))T ym . (22) 
t-voo z at It 



RATES 



The 2S/W model has six parameters: the four rates, 
7 and r , and the two positions x L / R . The inter- well 
rates we choose as twice the Kramer's rates 



V ± =2r K± , 



(23) 



where tk± are the Kramer's rates for passage out of the 
well to the right (+) and left (—)[!, Ha]. Explicit forms 
are given in Eq. (|A7|) . The factor two arises because, 
if the particle is localised in the right half of the well, 
the average time it will take to hit the right edge will be 
half that as if the particle was distributed over the whole 
well. The remaining parameters we fix by comparing the 
initial behaviour with that of the SE, Eq. ([1]). 

If we start the SE equation with the particle localised 
at a point x a and let the system evolve for a short time, 
then, to linear order in time t, the mean position of the 
particle evolves as a deterministic particle subject to a 
constant force f(x a ) = —u'(x a ). In contrast, the second 
cumulant evolves as if the particle experiences free diffu- 
sion. Evolution from x^m therefore gives the first two 
cumulants 

(x)~x L/R + f(x L/R )°trt; (x 2 ) c ~2D t. (24) 
Equating these results with Eq. (fT5)) and rearranging 



gives 



1 



2D Q 



7 = 7 = = / ^2 
2 (x R -x L y 

(x R - x L )f(x L ) = 2kT 

(x R - x L )f(x R ) = -2kT. 



(25) 



Equations (|25|) have to be solved numerically for x l and 
x R , which then give the rates 7 ± . 

With these choices, since 7 + = 7~, we obtain for the 
plateau values of Eq. (fT7| the values 

(x)^ = \{x L + x R ) ■ (z 2 )^ = \ (x L - x R f ; 

<* 3 }P lat =0; (x 4 )P' at = -I ( XL - x R ) 4 . (26) 
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FIG. 2. Comparison of diffusion coefficient D — 

limt->oo \-ji(x 2 {t)) c of Eq. (f22j) with the exact result of 
Ref. 0. The results of the 2S/W model agree well with 
the exact ones in the regime of deep valleys, i.e. uo <C kT 
and F <C F c , with F c the force for which potential valleys 
disappear. 



Furthermore, in the limit r ± /7 ;g> 1, the asymptotic cu- 
mulants read: 



-i) fc r- t, 



(27) 



which are identical with those of the one-state-per-well 
model (see Appendix |Aj 



VI. RESULTS 

Before discussing the short-time behaviour of the sys- 
tem, we begin by demonstrating that our 2S /W model re- 
covers the known long-time behaviour. In Fig. [2] we com- 
pare the diffusion coefficient D from the 2S/W model, 
Eq. ([22)) . with the exact result of Reimann et al. [Icj |. 
In accordance with expectations, our model works well 
provided that the valleys of the potential are sufficiently 
deep, Ait <C kT. This condition can be satisfied by en- 
suring that uo <C kT and that the applied force is less 
that the value of the force F c = uqtt / 'a at which the min- 
ima of the potential disappear. 

Turning now to the short-time dynamics, we consider 
first the second cumulant (x 2 (t)) c as a function of time in 
Fig. |31 where time is measured in units of r = Three 
distinct behaviours are seen. At the shortest times, the 
behaviour is essentially free diffusion with the second cu- 
mulant growing linearly. This is the behaviour which 
was fitted in our 2S/W model to obtain Eqs. (|25|) . As 
the particle begins to feel the influence of the wells, a 
plateau develops in the mean squared displacement. The 
duration of this plateau increases the higher the potential 
barrier between wells. This behaviour is well reproduced 
by the 2S/W model not just qualitatively, but quanti- 
tatively. We note that this behaviour was not explicitly 



10° 

io- 1 
io- 2 

IO" 3 



10 



-4 



.-2S/W,:: 
Un = 5kl 


SE 

n 


/ 






u = 
. u = 
- u = 


-- mm 

I5kl 
2Qkl 


n 
•i 
•n 
















- 











































1Q- 4 1Q- 3 IO' 2 



10 3 



t/r 



FIG. 3. The mean squared displacement (x 2 {t)) c for different 
amplitudes uo with no applied force (F = 0). The lines show 
the result from the 2S /W model, and the symbols those from 
numerical integration of the SE. For the parameters shown 
here, the 2S/W model reproduces the essential features of the 
short-time dynamics. The agreement with the numerically 
exact results increases with increasing well-depth (increasing 
M ). 



put into the model; rather, it emerges as a prediction of 
it. Fig. E] shows the height of the plateau as a function of 
the applied force F. The 2S/W approach works better 
the more pronounced the potential barriers, but even for 
uq = 5kT, the 2S/W theory and numerics agree quite 
well. Fig. |4] shows two different results for the plateau 
height extracted from the 2S/W model; one based on 
the full solution and one given by the approximate value 
of Eq. ([26]) . The coincidence of these results is a good in- 
dication that the mean-square displacement undergoes a 
genuine plateau, rather than something more like an in- 
flexion point, as is the case for e. g. uo/kT = 5 in Fig. [3] 
Not shown here, but the first cumulant (x) c shows simi- 
lar plateau behaviour, the height and duration of which 
is also well reproduced by the 2S/W model. 

At longer times, the particle breaks out of the confine- 
ment of its initial well, and starts diffusing through the 
potential. The second cumulant then starts increasing 
linearly again, this time with the diffusion constant of 
Eq. (f22|) and Fig. [2] The time at which the plateau goes 
over into diffusive behaviour is the cross-over time. Dur- 
ing the plateau and through the cross-over, the probabil- 
ity density slowly reshapes. This can be seen in Fig. [5l 
which shows the kurtosis, Eq. ((4]), as a function of time. 
By comparing Fig.[3]and Fig.[5]one can see that the time 
of the maximum deviation from Gaussian e.g. the maxi- 
mum of n(t) can be used as an approximation to the end 
of the plateau, as proposed in Ref. [T3|. It is remarkable 
that the dependence of this crossover time on uq and F, 
shown in Fig. [51 is given in common quantitatively. 
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FIG. 4. Plateau height (x 2 }^ 1Sjt as a function of driving 
force, F, calculated in three different ways. We have de- 
fined the plateau height as the value of (x 2 ) c at the time 
for which d(ln{x 2 ) c )/d(lnt) reaches its minimum value. The 
solid curves show the plateau heights obtained from the 2S/W 
theory, the symbols from integrating the SE and the dashed 
lines show the approximation of Eq. (|26[) . The 2S/W method 
underestimates the plateau heights slightly, becoming better 
as ua/kT becomes larger. Only when the two 2S/W results 
agree with one another does the cumulant really show a well- 
defined plateau. 



VII. CONCLUSIONS 

In this paper we have investigated the short-time dy- 
namics of a particle driven through a modulated poten- 
tial using a Markovian master equation approach based 
on a discretization of the spatial axis into two states per 
well. The remaining input into the theory are the transi- 
tion rates, which we have chosen by using Kramer's rates 
(inter-well transitions) and by comparison with the true 
short-time dynamics from the continuous (SE) approach 
(intra- well transitions) . 

A major advantage of our discretized two-state model 
is that relevant transport quantities, such as the mean- 
squared displacement, can be calculated analytically. In 
particular, we have derived explicit expressions for the 
heights of the plateaus in the mean-squared displace- 
ment at intermediate times, i.e., within the subdiffusive 
regime, and for the kurtosis, which signals the crossover 
to the diffusive long-time regime. The results from the 
two-state model are in excellent agreement with those 
obtained from numerical solution of the SE with and 
without external drive, provided that Au/kT 3> 1 and 
F < F c . This latter requirement is consistent with 
the expectation that (only) in this limit, the valleys in 
our continuous model can be approximated by almost- 
isolated wells. In the same range of parameters, the dif- 
fusion coefficient D (derived from the second cumulant), 
is found to be fully consistent with previous analytical 
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FIG. 5. The kurtosis n(t) of Eq. Q as a function of time 
for the same parameters as in Fig. [3] The dominant feature 
is a peak which occurs around the cross-over time. As be- 
fore, good agreement between 2S/W model and numerics is 
obtained for uo 3> kT. 



results of Rcimann et al. [Toj . It is precisely the high- 
barrier case, which is relevant in many experimental real- 
izations of colloids in modulated potentials (see e.g. 0]), 
as well as in the more general context of modeling sys- 
tems exhibiting cage effects (such as undercoolcd liquids 
[3(il ] and dense colloidal suspensions [13] )• For such sys- 
tems our "minimal" model for short-time dynamics could 
be particularly useful. 

We close this paper with some remarks on possible ex- 
tensions of our approach and its relation to other models. 
To start with, we have found that discrepancies between 
the master equation and the SE approach arise as Au/kT 
decreases. In this case, quantitative agreement can be im- 
proved by modifying the manner in which the rates are 
calculated. Firstly, expanding the full generating func- 
tion to linear order in time yields 

(x(t)) ~ x L/R ± 7 ± (:r R - x L )t T T T (a + x L - x R )t 
{x 2 (t)) c ~ y ± (x R - x L ft + r =F (a + x L - x R ) 2 t. (28) 



These values can be equated to Eq. (|24[) to match bet- 
ter with the SE. Secondly, the inter- well rates, Eq. ([2"5T) . 
can be improved by using higher-order corrections to 
Kramer's rate [{J. These latter corrections are partic- 
ularly important when the periodic contribution to the 
potential becomes steeper than the sinusoid considered 
here. 

A further point concerns the performance of our mas- 
ter equation approach for the prediction of higher-order 
cumulants. Whereas the complete behaviour of the first 
two cumulants is well reproduced by our 2S/W model, 
the agreement for the higher-order cumulants is not as 
good, and although the broad features of these k > 2 are 
reproduced, some fine details at very short times are not 
captured. To improve these higher cumulants it seems to 
be promising to augment the model with more states per 
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FIG. 6. The time of the kurtosis peak as a function of driving 
force F for different values of uq. Again, the agreement is 
good for deep valleys. 



well, an extension which is conceptually straightforward. 

Finally, we note that there is an connection between 
our model and other models describing subdiffusion. By 
tracing over the internal degrees of freedom, a nonMarko- 
vian master equation for the well-index n alone can be 
derived, involving a memory kernel which differs from 
a delta function (characterizing the Markovian case). 
Specifically, in the 2S/W model, the kernel in Laplace 
space has a single simple pole at — r, which describes the 
rate at which the intrawell degrees of freedom relax. It 
is well established [13] that such nonMarkovian kernels 
are intimately related [HI, [3i| to the waiting-time distri- 
bution in continuous-time random walks [29|, such that 
connexion can also be made with these latter methods. 
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particle in well n at time t. The position of the particles 
is given as x = na. The rate equation looks like 

= -{W + + W-)p {n) + W + p {n - 1] + W- P (n+1} (Al) 

where W ± are rates. We define the MGF and CGF via 



= g( x ,t) = J2 einx P (n) > 



such that the kth derivative 



^-pF( X ,i)|o = <« R W>c 



(A2) 



(A3) 



is the fcth cumulant of the well number. The correspond- 
ing position cumulants are simply: (x k (t)) c = a k (n k {t)) c . 
Fourier Transform of Eq. (|A1|) gives 



Q(x,t) = W( x )G(x,t) 



(A4) 



with W(x) = J2± W ± (e ±lx — 1). Solving this equation 
assuming that we start localised in the n = well, gives 
the CGF 

Hx,t) = W(x)t. (A5) 

Differentiating, we obtain the cumulants 

(x k (t)) c = a k [W+ + {-l) k W-]t. (A6) 

The cumulants of this model are therefore linear in time 
for all time, and there can be no caging effects in this 
one-state-per-well model. All odd cumulants are equal. 
All even cumulants are equal. 

Good agreement with the asymptotic solutions fioj is 
obtained in the regime Au/kT ^> 1 with the rates chosen 
according to Kramers rule, W ± = tk± with 0, |3(| 



r K-\ 
TK- 



£fcy/\u"{x mia )\\u»{x c )\ e -(«(«.)-«(*-*0)/w 
'^y/\W'(x mia )\\u"{x a )\ e -W^)-«(^))/p, 7 ) 



2-nkT 



where x m ; n , x a , x c are the positions of a potential mini- 
mum and the maxima to its left and right, respectively. 



Appendix B: CGF 



Appendix A: One-state-per-well model 

The simplest approach is to split the potential into 
a number of wells (labelled n = . . . , —2, —1, 0,1,2,... ) 
defined between subsequent maxima, and write down a 
rate equation for the probability p^ n >(t) of finding the 



We demonstrate here that Eq. (JTUJ) gives the position 
cumulants as stated in Eq. (fT2j). Consider the moment 
generating function corresponding to Eq. (|10[) : 

g( X ;t) = = ((^ o |e^ A e w «V O) (0)>>. (Bl) 

Taking the fcth-derivative, we find 



d{ix) 



d{ix) h 



Ao|e 4xA |p(x;*)» 



d(i X y 



]T((0ol£ e 4x(A+ "V n) W>> 



£«MA + n) V"' W» = ]>> a + x ^ k PL } + ( na + x n) k P { R = <**>™ 



(B2) 
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which is the fcth moment of the spatial distribution. The standard relation between MGF and CGF completes the 
proof. 
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